Scattering of elastic waves by a fracture zone containing randomly distributed cracks     

Jun Kawahara  Teruo Yamashita 《Pure and Applied Geophysics》1992,139(1):121-144

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.  相似文献   

Using numerical modelling to evaluate tsunami hazard near the Kuril Island     

L. B. Chubarov  Yu. I. Shokin  K. V. Simonov 《Natural Hazards》1992,5(3):293-318

A sequence of computer experiments is used to study questions concerning the tsunami problem as a quantitative estimate of tsunami danger, detailed geographical tsunami classification, determination of the parameters of critical tsunami waves, and the conditions of their development. We call a wave critical, if its impact on the coast is most hazardous.Using the Middle Kuril Island as an example, we present the results of a computer experiment which includes determining the wavefields on the shelf and estimating the effects connected with the deep-water Bussol and Diana Straits.Numerical simulation of tsunami waves of different sources permits the assessment of the extent of tsunami danger in different areas of the coastal zone of Simushir Island, depending on the location of the focus zone and their geometry.The major singularities of the wavefield arise in the zones of the deep-water straits. The distribution of the amplification factors is determined by both the global parameters of the wavefields and the local properties of individual harbours. The results obtained for a particular harbour in the northern part of Simushir Island, formed the basis for the quantitative estimate of tsunami danger for this area.  相似文献   

Transmission and reflection of Rayleigh waves at a thin low velocity vertical layer: A laboratory and theoretical study     

P. F. Biagi  T. S. Gegechkory  E. N. Its  P. Manjgaladze  V. Sgrigna  D. Zilpimiani 《Pure and Applied Geophysics》1990,133(2):317-327

A three-dimensional ultrasonic modelling was used to study the Rayleigh waves reflected from a thin low velocity vertical layer and those transmitted through the same zone. The reflection and transmission coefficients of the Rayleigh waves have been calculated by an approximate method based on the Green function technique. Nonrigid contact boundary conditions were used for the low velocity layer modelling. A rather good correspondence between the theoretical and experimental absolute values of the reflection and the transmission coefficients was obtained. Some discrepancies between experimental and theoretical results, in particular with regard to phase shifts, may be explained by a certain inadequacy of the theoretical model to the low velocity layer investigated by the ultrasonic modelling.  相似文献   

弹性波的三维有限元模拟   总被引：5，自引：0，他引：5       下载免费PDF全文

邓玉琼  张之立 《地球物理学报》1990,33(1):44-53

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.  相似文献   

The sequence of sediment behaviour during wave-induced liquefaction   总被引：4，自引：0，他引：4  

B. MUTLU SUMER  FIGEN HATIPOGLU  JØRGEN FREDSØE  S. KAAN SUMER 《Sedimentology》2006,53(3):611-629

This paper presents the results of an experimental investigation of the complete sequence of sediment behaviour beneath progressive waves. The sediment was silty with d ₅₀ = 0.060 mm. Two kinds of measurements were carried out: pore-water pressure measurements (across the sediment depth), and water-surface elevation measurements. The process of liquefaction/compaction was videotaped from the side simultaneously with the pressure and water-surface elevation measurements. The video records were then analysed to measure: (i) the time development of the mudline, (ii) the time development of liquefaction and compaction fronts in the sediment and (iii) the characteristics of the orbital motion of the liquefied sediment including the motion of the interface between the water body and the sediment. The ranges of the various quantities in the tests were: wave height, H  = 9–17 cm, wave period, T  = 1.6 sec, water depth = 42 cm, and the Shields parameter = 0.34–0.59. The experiments reveal that, with the introduction of waves, excess pore pressure builds up, which is followed by liquefaction during which internal waves are experienced at the interface of the water body and the liquefied sediment, the sequence of processes known from a previous investigation. This sequence of processes is followed by dissipation of the accumulated excess pore pressure and compaction of the sediment which is followed by the formation of bed ripples. The present results regarding the dissipation and compaction appear to be in agreement with recent centrifuge wave-tank experiments. As for the final stage of the sequence of processes (formation of ripples), the ripple steepness (normalized with the angle of repose) for sediment with liquefaction history is found to be the same as that in sediment with no liquefaction history.  相似文献   

A laboratory experiment to monitor the contact state of a fault by transmission waves   总被引：1，自引：0，他引：1  

Naoto Yoshioka  Koji Iwasa 《Tectonophysics》2006,413(3-4):221-238

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.  相似文献   

Rock fracturing by explosive energy: review of state-of-the-art   总被引：1，自引：0，他引：1  

M. R. Saharan  H. S. Mitri  J. L. Jethwa 《Fragblast: International Journal for Blasting and Fragmentation》2006,10(1):61-81

A study of the dynamic rock fracture initiation and propagation due to explosive energy is presented through a detailed state-of-the-art review. Explosive energy dissipation in crushing and fracturing is examined and the various means to enhance the explosive energy utilization for dynamic rock fracturing are reviewed. The study highlights the need for a better understanding of the dynamic fracturing process particularly in the presence of in situ stresses in the rock mass.  相似文献   

The effective elastic thickness of the India Plate from receiver function imaging, gravity anomalies and thermomechanical modelling   总被引：2，自引：0，他引：2  

György Hetényi  Rodolphe Cattin  Jérôme Vergne  John L. Náb&#;lek 《Geophysical Journal International》2006,167(3):1106-1118

  相似文献   

S-wave velocity structure, mantle xenoliths and the upper mantle beneath the Kaapvaal craton   总被引：1，自引：0，他引：1  

Angela Marie Larson  J. Arthur Snoke  David E. James 《Geophysical Journal International》2006,167(1):171-186

  相似文献   

Re-evaluation of focal depths and source mechanisms of selected earthquakes in the Afar depression   总被引：1，自引：0，他引：1  

L. Hagos  H. Shomali  R. Roberts 《Geophysical Journal International》2006,167(1):297-308

We present a stepwise inversion procedure to assess the focal depth and model earthquake source complexity of seven moderate-sized earthquakes  (6.2 > M _w > 5.1)  that occurred in the Afar depression and the surrounding region. The Afar depression is a region of highly extended and intruded lithosphere, and zones of incipient seafloor spreading. A time-domain inversion of full moment tensor was performed to model direct P and SH waves of teleseismic data. Waveform inversion of the selected events estimated focal depths in the range of 17–22 km, deeper than previously published results. This suggests that the brittle–ductile transition zone beneath parts of the Afar depression extends more than 22 km. The effect of near-source velocity structure on the moment tensor elements was also investigated and was found to respond little to the models considered. Synthetic tests indicate that the size of the estimated, non-physical, non-isotropic source component is rather sensitive to incorrect depth estimation. The dominant double couple part of the moment tensor solutions for most of the events indicates that their occurrence is mainly due to shearing. Parameters associated with source directivity (rupture velocity and azimuth) were also investigated. Re-evaluation of the analysed events shows predominantly normal faulting consistent with the relative plate motions in the region.  相似文献